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Cluster variables have recently revolutionized numerical work in certain models involving fermionic variables. 
This novel representation of fermionic partition functions is continuing to find new applications. After describing 
results from a study of a two dimensional Hubbard type model that confirm a superconducting transition in the 
Kosterlitz-Thouless universality class, we show how a cluster type algorithm can be devised to study the chiral 
limit of strongly coupled lattice gauge theories with staggered fermions. 



1. INTRODUCTION 

During the last few years a new class of 
fermion algorithms have emerged. The essential 
progress is a result of our ability to rewrite cer- 
tain fermionic partition functions as a sum over 
configurations of bond variables with positive def- 
inite weights [0J|], i.e., 

Z = Y, W[b] SiiH[6] (1) 

[b] 

where W\b] > is the weight of the bond configu- 
ration \b] and Sign [6] > is an entropy factor that 
takes into account degrees of freedom other than 
the bond variables. Typically, the Pauli principle 
is encoded in the topology of clusters formed by 
lattice sites connected through the bonds. Clus- 
ters also carry a variety of interesting physical in- 
formation. For example, sizes of certain clusters 
are related to condensates, the squares of the sizes 
of clusters yield susceptibihties. Further, clusters 
are useful in building efficient algorithms close to 
critical points where the correlation lengths di- 
verge since they allow non-local updates with a 
reasonable acceptance. This property has helped 
in studying critical phenomena in fermionic mod- 
els with unmatched precision. 

2. SUPERCONDUCTIVITY 

The recent success of cluster methods in 
fermionic systems originates from the Hamilto- 
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nian approach. One starts by writing the parti- 
tion function in terms of configurations of fermion 
occupation numbers whose Boltzmann weight is 
not guaranteed to be positive definite. By intro- 
ducing additional bond variables which help in 
the integration over the fermion occupation num- 
bers, the partition function of certain models can 
be rewritten in the form of eq. (|l]). Configura- 
tions with clusters of a specific topology, called 
merons, then yield Sign[6] = 0. 

Recently, superconductivity in a two dimen- 
sional attractive Hubbard type model was studied 
using the meron cluster approach. The fermion 
pairing susceptibility (x) is a useful observable. 
It is expected to satisfy the finite size scaling for- 
mula 

^^)-{ Const. T>T: 

if the superconducting transition belongs to the 
Kosterlitz-Thouless (KT) universality class, with 
< ri{T) < 0.25, 7?(rc) = 0.25 and ?7(0) = 0. In 
the specific model studied, (x) turns out to be a 
sum over the square of the size of each cluster in 
the zero meron sector and the product of the size 
of the merons in the two meron sector. Figure |^ 
shows a plot of the susceptibility as a function of 
lattice size for various temperatures. Consistency 
with KT predictions is clear. 

Another observable relevant to the study of 
superconductivity is the spatial fermion winding 
number susceptibility (W^) . Although this is dif- 
ficult to evaluate with conventional algorithms, it 
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strained. Local metropolis updates which can be 
formulated in the massive case become exponen- 
tially inefficient in the chiral limit. Here we argue 
that cluster representations of the MDP systems 
yield useful algorithms directly in the chiral limit. 

To understand the cluster representation con- 
sider for simplicity the strongly coupled t/(l) 
gauge system. The partition function in this case 
is given by the number of closely-packed-dimer 
(CPD) configurations on a lattice. A typical CPD 
configuration in two dimensions is shown in Fig. 



Figure 1. Pairing susceptibility as a function of 
system size for various temperatures. 



is relatively straight forward in the meron cluster 
approach. Each cluster can be assigned a spa- 
tial fermion winding number. The susceptibility 
then turns out to be the sum over the square of 
each cluster's spatial winding number in the zero 
meron sector and the product of the spatial wind- 
ing number of meron clusters in the two meron 
sector. In the infinite volume limit below the crit- 
ical temperature, one can combine known results 
to obtain 2'K'q{T){W^) = 1. Results again show 
consistency with this expectation. Preliminary 
results from this study was presented in |^ and 
the final analysis in Q. 

3. CHIRAL SYMMETRY BREAKING 

Although the recent success has been applied to 
Hamiltonian models of chiral symmetry breaking 
Q, cluster methods are applicable to more con- 
ventional Lagrangian models as well. For exam- 
ple, consider strongly coupled lattice gauge the- 
ory with massless staggered fermions in which the 
[/(I) chiral symmetry is expected to be broken 
spontaneously in four dimensions j^. This re- 
sult was obtained by mapping the massive model 
into a statistical mechanics of monomer- dimer- 
polymer (MDP) system with positive definite 
Boltzmann weights and extrapolating the results 
to the chiral limit. Unfortunately, as far as we 
know, it has been difficult to devise algorithms in 
the chiral limit where the systems become con- 



Figure 2. A two dimensional CPD configuration 

Such configurations are also of interest in sta- 
tistical mechanics and play an important role in 
the solution to the 2-d Ising model Q]. The chi- 
ral symmetry of staggered fermions is manifest 
in this representation by the fact that the chi- 
ral condensate vanishes since it is impossible to 
find a CPD configuration with one defect (one 
site has no dimer lines attached to it). The chiral 
susceptibility on the other hand is non-zero and 
proportional to the ratio of the total number of 
CPD configurations with two defects (two sites 
are not connected by dimers) and the partition 
function. 

It is possible to extend CPD configurations 
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Figure 3. Rules for extending the CPD configu- 
rations to include additional bond variables. 



susceptibility with system size. Surprisingly, al- 
though the data is not inconsistent with the pres- 
ence of long range correlations, the susceptibility 
does not seem to follow the predictions of eq. (||) . 
This puzzle is currently being investigated along 
with extensions to higher dimensions. 
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to configurations of loops made up of bonds 
which include the original or "filled" dimers (rep- 
resented here by "solid" bonds) and "empty" 
dimers (represented by "dashed" bonds) such 
that the partition function can be expressed as a 
sum over weights of new loop configurations. Fig- 
ure 1^ shows the rules of one such extension in two 
dimensions. Each shaded plaquette of the CPD 
configuration of Fig. ^ carries one of the seven 
plaquette configurations given on the left side of 
equations in Fig. |3[ It is easy to check that all 
constraints are satisfied if each loop is made up of 
a repeating sequence of filled and empty dimers. 
The usefulness of the loop variable is that a dimer 
system can be updated by "flipping" a loop where 
filled dimers are emptied and vice versa. The ac- 
ceptance of such a flip is reasonable and leads to 
a useful algorithm. 

The chiral susceptibility gets contributions 
when a part of the loop is flipped and can be 
measured easily along with the update. The al- 
gorithm was first applied to the two dimensional 
model. Although a t/(l) chiral symmetry can- 
not spontaneously break in two dimensions, long 
range correlations can arise as predicted by the 
Kosterlitz-Thouless universality class as discussed 
in the previous section. Figure ^ plots the chiral 



Figure 4. Chiral susceptibility as a function of 
lattice size in two dimensions. 
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